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Q-i Abstract. The boundary homogenization for Dirichlet problem 



c/3 



(system), is the subject of study in this paper. Using methods 
from Fourier analysis we prove LP convergence for solutions to 
the system 

div(A(x)V« f (j)) = 0, inD, u E (x) = g(x,x/e) ondD, 

to a limit function Uo, for bounded, uniformly convex, regular 
domains D c ]R d (d > 2), and regular operators and boundary 
data. The outcome of our paper is the following convergence rate 
for all 1 < p < oo 

ie 1 ' 2 ?, d = 2, 

Hw £ - «olb(D) <C P (ellnel) 1 ^, d = 3, 
\e ll P, d > 4, 

which we prove is (generically) sharp for d > 4. Here Mo is the solu- 
tion to the averaging problem. As a consequence of the proof of the 
convergence result, we also obtain estimates for the convergence 
rate for homogenization of Neumann problem with oscillating 
boundary data. 



1. Introduction and main result 

Recently, the current authors obtained pointwise convergence rate 
for solutions to the boundary layer, and homogenization problem for 
Dirichlet system, with rapidly oscillating data. From our pointwise 
estimate one could easily obtain U convergence rate for solutions, of 
the oscillating problem to the effective limit. That the effective limit 
exists, is a consequence of the assumption of convexity and regularity 
of the domain, as well as the smoothness of the operator and the 
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boundary data (see [LSI, [ASSj ). Nevertheless, the results in flASS] , 



not being optimal, raise naturally the question of finding the optimal 
rate for U convergence. Similar attempt has been made in several 
papers, where sometimes strong assumptions have to be made in 
order to obtain some reasonable results. The latest results for the 
convergence rate can be found in D. Gerard- Varet and N. Masmoudi 



HGML where they also consider oscillating operators. This in its 
own turn implies substantial difficulty as the interplay between the 
oscillating operator and the data gives rise to new unknown terms. 

In this paper we shall strengthen our results on V convergence 
rate in HASS1 , which we prove to be optimal in R rf for d > 4, and close 
to optimal in R 3 . A partial result concerning optimality in dimension 
two is also considered. Our technique uses Fourier analysis meth- 
ods and depends heavily on the regularity of the operator, and the 
boundary data, as well as on the regularity, and uniform convexity 
of the domain. The method is straightforward and computational 
analysis, using estimates of the Poisson kernel, and is based on the 
study of oscillatory integrals with singular kernels. 

Let D be a bounded and uniformly convex domain in R rf (d > 2), 
and T be its boundary. Let also A{y) = (AJ (y)), 1 < a, ft < d, 1 < 

i, j < N be an R N2x<i2 -valued function defined on R rf , and g be an (Un- 
valued function defined on D x TP 1 . We study asymptotic behavior 
of solutions to the following problem: 

(£u e (x) = 0, in D, 

\u E (x) = g(x,x/e), onT, 

where e > is a small parameter, and using the summation conven- 
tion of repeated indices the operator L is defined as 



(1) 



<r ^ d 



dx a 



y dxP 



-div [A(x)Vu] , 



where u = (« l7 u N ) and 1 < i < N. 

For (d) we consider the corresponding homogenized problem 



Lu (x) = 0, in D, 
Uo{x) = g(x), on T, 



(2) | 

where g(x) = f g(x, y)dy. 

It is interesting to note that our methods developed for homoge- 
nization of the Dirichlet boundary value problem, enable us to treat 
the case of a similar problem with Neumann boundary value. For 
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the operator £ defined above, and for some function F e : M. d — > C N 
consider the following problem 



(3) 



Lu £ (x) = F e (x), in D, 



^(x) = g(x,x/e), onT, 



where (x) = n a {x)A a J.{x)^, 1 < i < N, denotes the conormal 

derivative, and n(x) is the outward unit normal to T at the point x. 
Here for each e > one chooses F e so that the compatibility condition 

J F £ (x)dx = J g(y, yle)do{y) holds true, and sup ||F e |L < oo. 

d r e>o 

Remark 1.1. It is noteworthy that our results can be formulated in terms of 
double and single layer potentials, regardless of solutions to partial differen- 
tial equations. In this regard the method presented here gives the advantage 
of being applied to problems in potential theory, and double/single layer 
potentials with oscillating densities. 

1.1. Assumptions. We make the following assumptions: 

i (Periodicity) The boundary vector-valued function g is 1- 
periodic in its second variable: 

g(x, y + h) = g(x, y), Mx e D, My e K d , V/z e X d . 

ii (Ellipticity) There exists a constant c > such that 

iii (Convexity) We assume that T is uniformly convex hypersur- 
face, that is all its principal curvatures are bounded away from 
0. 

iv (Smoothness) We suppose that the boundary value g in both 
variables, function F e in problem ^ for each e > 0, all ele- 
ments of A and the surface T are sufficiently smooth. Here we 
do not aim to obtain the optimal smoothness, but rather focus 
on the method itself. 

v (Symmetry) In case of problem © we assume A = A*, i.e. 



The main results of this paper are the following. 



a a ?(y) = aj(y) for all 1 < i,j < N and 1 < a,$ < d. 



Theorem 1.2. (Dirichlet problem) Let u e be a solution to the system ® 
and Uq to that of © under assumptions (i)-(iv). Then, for all I < p < oo 
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one has 



u \\ 



LP(D) 



e 1/2 P, d = 2, 

(ellnel) 1 ^, d = 3, 
eVp, d>4. 



Our efforts to get rid of the log-term in the case d = 3 did not 
give any result, partly due to the integrals containing the highest 
derivatives of the kernel (see e.g. integral A 6 in the proof of the main 
theorem). 

Next, we consider the question of optimality of the U convergence 
rate provided by Theorem II .21 In particular we prove that in dimen- 
sions greater than 3 the convergence rate obtained in Theorem II .21 is 
sharp. For simplicity we will consider the case of simple equations 
rather than system. 

Theorem 1.3. (Optimality) Let N = 1, and D, T, and L be as before. Let 
u e be a solution to the system © and w to that o/© under assumptions 
(i)-(iv). Then for each 1 < p < oo there exists a constant C p independent of 
e, such that 

\\U £ - U \\ L P(D) > C p £ Vp \\g - glli-fi*). 



Theorems II .21lL3l imply that the convergence rate of homogeniza- 
tion of the Dirichlet problem is optimal when d > 4. 

As an application of the proof of convergence result for the Dirich- 
let problem, we consider homogenization of the Neumann problem, 
with oscillating boundary data. 

Theorem 1.4. (Neumann Problem) Let u e be a solution to the system © 
under conditions (i)-(v), and u be a solution to the same problem where 
the boundary value g is replaced by g, and F £ is replaced by some smooth 
function F to fulfill the compatibility condition. Denote 

v e (x) = u E (x) - jij- J u £ (y)do(y)- J N(x,y)F e {y)dy, 
r d 

and let Vq be the term corresponding to homogenized problem. Then for 
d>3 and 1 < p < oo one has 



Vo\\ 



U'(D) 



<c v 



e 1 ^, d = 2, 

e y P, d = 3, 

e 3 ' 2 ?, d = 4, 

e 2/ P | lne |i/ P/ d > 5 _ 
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The reader may wonder about the behavior for the Neumann prob- 
lem, versus Dirichlet above, and the logarithmic term that appears 
for all dimensions greater than 4 in the Neumann case! This non- 
optimality with log-term, appears exactly in our computation of the 
integral which includes the highest singularity of the representation 
kernel (see e.g. integral As in the proof of the main theorem). 

The following is an example of problem |3l for which the conver- 
gence rate of its solutions is determined by its boundary data. 

Example 1.5. For each e > take F E = p- J" g(y,y/e)do(y), and F = 

r 

jSf f ~g(y)do(y)- Since g is sufficiently smooth function, and T is a smooth 
r 

and uniformly convex hypersurface, after expanding g into its Fourier series 
with respect to the periodic variable, and applying the principle of stationary 
phase (see (Sj, chapter VIII, Theorem 1) on each summand we have 

\F E - F | < Ce id ~ 1)/2 . 
Using this, and Lemma \T7W\ we obtain 

J N(x, y)(F £ - F )dy < Ce {d ~ 1)/2 J \N(x, y)\dy < Ce (rf " 1)/2 , 

D D 

where C is independent ofxeD, and e > 0. Combining this last estimate 
with Theorem FT4\ for each 1 < p < oo we arrive at 



u ° ~ W\ I ^ £ ~ u °" >do ^ LP(D) - Cp 



'e 1 ^, d = 2, 

e y P, d = 3, 

e 3 ^, d = 4, 

The example shows, that we will have the same picture, if we take some 
smooth and one periodic function F(x), and proceed by taking F £ (x) = 
F(x/£), andF = J F{x)dx. 

Theorem 1.6. (Gradient of Neumann solutions) Let v E and v be as in 
Theorem H~?l Then for 1 < p < oo, and any k > where x < l/2p for 
d = 2, and k < 1/pfor d > 3, one has 

\N(v e - v )\\ mD) < C p , K £ K . 

The authors believe the result of Theorem [L2] is optimal in dimen- 
sions two and three as well. The optimal rate also indicates that one 
also expects that the problem with both oscillating coefficients and 
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boundary data should have a similar rate, at its best. An evidence 
for this can be observed in IIKLS1I1 where using the so-called Dirichlet 
correctors the problem with both oscillating operator and boundary 
data is transformed to a new problem with only oscillating Dirich- 
let data, where the difference between the solution of the original 



problem and the new one (see [KLS1J, Theorem 3.9) is similar to the 



estimate in our Theorem 11.21 Although this transforms the general 
problem to one with fixed operator, we can not apply our method 
directly to homogenize the new boundary data, since the new oscil- 
lation term is not explicit, and involves the Dirichlet correctors. 

We also note that our method depends on the regularity of the 
Poisson kernel and its derivative. The Poisson kernel has the rep- 
resentation P(x, y) = n(y)A(y)V y G(x, y), with n{y) the outward unit 
normal to T at y, G the Greens function and A the matrix represent- 
ing the operator. This implies that oscillating operator P e (x,y) = 
n(y)A(y/e)V y G e (x, y), has derivatives with the factor 1/e in front of it 
which causes our method to fail in the case of oscillating operators. In 
the paper of IIGMI1 this is taken care of by approximating the surface 
by hyperplanes, scaling techniques and the use of correctors. These 
steps of IIGM1 naturally cause the loss of the optimality in speed. 
Finally we refer to some standard texts in PDE IADN1I , BADN2I1 , 



homogenization theory |AA| , |BLP] , |JKO| , and oscillatory integrals 
llSl for the background. 

1.2. Preliminaries. Before proceeding to the proofs of the main re- 
sults, we need the following lemmas. 

Lemma 1.7. ( HASSL Lemma 2.1) Let P(x,y), where x e D and y £ T, 
be the Poisson kernel for the problem ((T]) under assumptions (i)-(iv). Then 
for each a = {a\, ...,ai) e Z+ there exists a constant C a depending on a, 
domain D, and operator £, such that 

\D«P(x,y)\<C a ] ^-±^ , xeD,yeT, 



|P(x, y)| < C 0T:r ^ R , xeD,yeI, 



d(x) 
\x - y\ a 

where \a\ = |«i| + ... + \aA, and d(x) is the distance of the point xfrom the 
boundary T. 



The first estimate follows from HASSH and the second (with dis- 
tance) from IIKLSIB equation (3.1). 



Using this Lemma we can establish uniform bounds with respect 
to x G D of the surface integral of the modulus of P(x, y), which we 
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will use later on. Fix x e D, and take r > 0. Due to convexity of the 
domain the level set E r := \y e T : \x — y\ = r} is contained in a base of 
a circular cone with vertex x and side length r, and hence there exists 
a constant C depending on the domain D so that "Kd-iiL,) < Cr d ~ 2 . 
Using this observation and applying the second estimate of Lemma 
ll.7l we obtain 

diam(D) 

(4) J \P(x,y)\do(y) < Cd(x) J < Cd(x) J T -^dr < C, 

r r d(x) 

where the constant C is independent of x e D. 

Lemma 1.8. ( lASSl , Lemma 2.3) Iff e C k (T d ) and jS e R, f/ien 

provided k + $ > d/2, where c m (f) is the m-th coefficient of the Fourier 
expansion of f, \\m\\ denotes the Euclidean norm of the lattice point m, and 
\a\ = ot\ + ... + ad- 

An immediate consequence of Lemma lL8l is the following result. 

Lemma 1.9. Let ft eK,Dbea compact subset o/R rf , and f a given function 
f(x, y) : DxT rf -> R, which is periodic in its second variable. Suppose that 
for all a e Z+ with \a\ < k, Dyf(x, y) exists and is continuous on Dx T d . 
Then 

\c m (f;x)\ 



E 



IMP 

m*0, meZ d 



< C k+ p,f, VX € D, 



provided k + f5 > d/2, where c m (f;x) is the m-th coefficient of the Fourier 
expansion of f at the point x. 

By N(x, y) denote the matrix of Neumann functions for operator £ 
in D (see IIKLS3I1 for definition). We prove the following lemma. 

Lemma 1.10. Under the assumptions (i)-(v) and d > 3, for each a = 
(cti, a,]) s there exists a constant C a such that for allx e D and y eT 
one has 

1_ 

\x - y\ d+ ^ a 

and 

\D a y V x N{x, y)\ < C a - \-r-r, (II), 



\D" y N(x,y)\<C a - (I), 



\ x _ y\d+\a\-l' 
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where \a\ = \a\\ + ... + \atd\- 

For d = 2 one obtains the estimates of similar type where d + \a\ — 2 is 
replaced by \a\ + x, in the first term, and d + \a\ — 1 is replaced byl + \a\ + k, 
in the second term. Here k can be take to be any positive constant. 

The estimate for Neumann function in R 2 is not optimal in the 
literatures; at least we could not find a reference. One should expect 
the logarithmic estimate, but the best estimate we found is that in the 
above lemma (see Remark 5.7 in j|KLS3fl ). 

Proof. The case when \a\ < 1, under weaker conditions on operator 
and domain was treated in HKLS3H . The case of \a\ = 2, or even higher 
orders can be done by a scaling argument along with up to boundary 
uniform regularity for solutions to Neumann problems; see Lemma 
2.1 in IIASSI1 for a similar treatment for the Poisson kernel. □ 

An easy consequence of this lemma is the following bound on the 
gradients of u e . 

Lemma 1.11. Let u e be a solution to the problem ©. Then for each x > 
there exists a constant C K independent of e such that 

|Vw £ (*)l < C K -^- f Vx G D. 
d K {x) 

Proof. The following representation is known 

(5) u cW-^\§ u E (x)dx = J N(x,y)F £ (y)dy + J N(x,y)g £ (y)do(y). 



J \v x N(*,y)\dy + J \v x N(x,y)\do(y) 



Using the uniform boundedness of F e and g £ with respect to e > 0, 
from (0 we obtain 

(6) |Vtt e (*)| < C 

The volume integral in ^ is bounded by virtue of estimate (II) in 
Lemma 11.101 For the surface integral, again the estimate (II) of 
Lemma [1.10[ gives us 

do(y) C f do{y) 



Jlv MxM) <cJJ^.<S- j f 



x-y\ 



d-l-K ' 



The last integral is uniformly bounded with respect to x by some 
constant depending on k, hence we obtain the result. □ 
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2. Proofs of convergence results 

Proof of Theorem 11.21 We divide the proof into some steps. 

Step 1. Reduction to local graphs. Let z e T and r > be small. Then 

there exists an orthogonal transformation Ji such that 

(7) (fl(r - z)) n B(0,r) = {(y'^(y')) : ly'l < 10r} n B(0,r), 

where y' = (yi, y rf _i), i/>(0) = V</>(0) = and ^jf(0) = a jf ] = 

1,2, ...,d - 1, with > 0, and ^(0) = 0, for i * j. Also \D a \p\ < C a , 
and < c < a x < a 2 < ... < fl^-i < C. We also have 

(8) Kily1^IViKy')l< WL 

where and do not depend on z. Now choose 6 > so small that 

(a) 5 < Tj^j and K\b < 1, 

(b) ©holds for |y'| < J^S, 



(c) 



^ io5m for z ' * / and 



d 2 4> _ 



< ^ for |y'| < 1006 when 



j = \,2,...,d-\, 

(d) |n'| < K\b implies that there exists a unique |y'| < <5 so that 
V^(y') = n'. ' 

We remark that 5 is a constant that does not depend on z. We 

M 

have T c U B(z, ±L<5) where L = J^, hence r c U B(z fc , ^<5) for some 

zeT 2 k=l 

M 

z 1 , ...,z M G T. We take a partition of unity X <Pk = 1 on T, where 

k=l 

supp((p k ) c B(z ,c ,LS), and <p fc e C°°. Set B k = B{z k ,Lb). Recall that 
g(x) is the average of g on the unit torus with respect to its periodic 
variable, also denote g e (x) := g(x,x/e). 
We have 

M r M 

u e (x) - u Q {x) = I p (x, y)[geiy) - giy)\(pk{y)do{y) ■= ^ h, 

k=l y k=l 

where 

4 := J" P(x, y)[g £ (y) - g{y)\(p k {y)do{y). 

Step 2. Reduction to volume integrals. Set z = y - z k , then 

J* = J P(x, z k + z)[g £ - g](z k + z)cp k (z k + z)do(z). 

(TnB k )-z k 
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We have (r n B k ) -z k = (T- z k ) n B(0, Lb). By setting y = ftz we obtain 



h= J P(x,z k + 5C l y)[g £ - g](z k + 1C l y)<p k (z k + l 3C 1 y)do{y). 

%(T-z k )r\B(0,L5) 

By and (a) we may assume that 

tt(T - z k ) n B(0, Lb) = {(y', xl>{y')) : \y'\ < 1006} n B(0, Lb), 
and hence 

\y'\<L5 

cp k (z k + %-\y'^{y'))){l + \V^(y')\ 2 f 2 dy'. 

Step 3. Reduction to oscillatory integrals. Since g is one periodic in 
its second variable and sufficiently smooth, we have 

g(x,y)=Y j c m (x)e 2nimv , 

meZ rf 

and hence 

g E (x) = £ c m {x)^ x , 

meZ d 

where c m : R rf — > C N for each ra e Z rf . Using this and orthogonality 
of 31 we have 

(9) [g £ -g](z k + K-\y',W))) = 

£ c m (z* + ^- 1 (y^^(y')))e 2ra ^ I • 2 "^' 2 ™ l<:Rm ' (!/ ' W)) ^ 

where < •, • > denotes the scalar product. By setting n := Jim and 
n = \n\(n',nd) with \{n' ,n^)\ = 1 from ^ we obtain 

[g £ -g](z k + X- 1 (y',ip(y'))) = 

c,n(z k + %-\y', ^{y')))e lnil - m - zk e lniAH y'\ 

where 

F(y') = n' -y' + n^(y'), 
and A := y = since "R is orthogonal. Next, by setting 

o fc (y') = <p,(2 fc + x-V/ + Wy')l 2 ) 1/2 , 
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and 

\y'\<L5 

we obtain 

141 < Yj M- 

Step 4. Decay of 1^, We split the study of decay of the integrals \ m 
into two cases. 
Case 1. \n'\ > Krf/2. 

We have VF(y') = n' + n d Vxp(y'). Then |n'.| > K 1 5/2d for some 
1 < j < d — 1, hence by © on supp(Oj;) we have 



(10) 



dF 



> — - - K 2 \y | > - K 2 Lo = — : > 



2d 



2d 



2d 4d 



Ad 



Now integrating by parts in \ m in the ;-th coordinate twice, by virtue 
of ((10)) and Lemma HT7I for all x e D we conclude 



(11) \h,m(x)\ < CA" 



\y'\<L6 



\c m \ + 



3c„ 



r9l/ ( 



+ 



d 2 C 



\x-z k -%-i(y',xl>(y'))\ 



AMd-l+2 



Recall that d(x) is the distance of x e D from the boundary of D, and 
set 



(12) 

Now observe that 
(13) 



D f = ji G D : d(x) > e}. 



\x - z k - ^{y' ,^){y'))\ d+1 



dx 



r d-i c 

-r-rdr = —. 

r d+l £ 



D £ |h;|>£ 

Combining this and ((11) we obtain 



fi,,„ 



(x)\dx < CA 2 e 



d, 



-I 

\y'\<U> L 



dCn 



+ 



d 2 c n 



(z k +X-'(y',4>(y')))dy'. 



Now taking into account the smoothness properties of g and applying 
Lemma [L9l to g and to its derivatives to sum up c m and its derivatives, 
from the last estimate we obtain 

(14) Yj Wb(D 8 )<Ce. 

m*0,meZ d 
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Case 2. \n'\ < 

Since 5 > is small and \(n' , = 1 we have |n^| > 1/2 and hence 

Ki<5 



n' 



< 



Ik 



By (d) there exists a unique y' with |y'| < 5 and Vi/Xy') = — jj-. Clearly 
VF(y') = 0, and using (c) we arrive at 

d 2 ^ t 1 99 



*- 









From the latter it follows that 



-(y> + se ? 



99 

> fli s , s < 46, 

200 



where s is scalar and e'. is the /-th unit vector of R . By (c) for z j 
we have 



dF /3F\ 
dyi [dyjj 



rid 



d 2 ^ 



dy i dy ] 



< 



lOOOrf' 



from which we obtain 
(15) 



> c\zl 



for z' G Qj and y' + z' G supp(Oj;) where 

1 



e y = {z' g iR 1 



rf-i 



|z)l > — ;=|z% 



l,2,...,d-l. 



Clearly the cones C ; cover R rf_1 . For j = 1,2, d - 1 there exists &>; 
supported in Cy, smooth away from the origin and homogenous of 
degree such that 

d-l 

Y J OJ i (z') = h Vz'^0. 

7=1 

Now fix a nonnegative function h G C°(R d_1 ) such that h(y') = 
for |y'| > 2 and h(y') = 1 for |y'| < 1. Setting y' = y' + z' and 
z* := z k + ^ _1 (y' + z', i/Xy' + z')) we obtain 

= J c m (z)P{x,zy lniA ^' + ^ k {y' + z')dz'. 

\y'+z'\<U 
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Set 

It 



l{ m = J h{e- ll2 z')c m {z)P{x, z)e lmK ^>'^<$> k {y> + z')dz' , 

\y'+z'\<L8 



and 



7 2 

k,m 



j (1 - h{e- ll2 z'))c m {z)P{x, z)(? niAF{y ' +z,) ® k {y' + z')dz' , 

\y'+z'\<L5 

so that J; C/ ,„ = 1} + J? . It follows from Lemma [L7l that 

dx 



x-y\ 



uniformly for y G T and e > 0, which together with the smoothness 
condition on g and Lemma [L9l gives 

X / KM** ^ c L / M^kfa' * c^~ 1)/2 - 

For the second part we have 7? = ]T I, where 

-i fc,m *— 1 k,m 

7 S, = J a; / (z')(l-^(e- 1/2 z / ))c m (z*)P(x,z*)e 2ra ' AF 0,(^+ Z ')rfz'. 
Now integrating by parts with respect to Zj in it' twice we obtain 



(16) \il(x)\<CA- 2 - 



f h t {i^ [ i^ (z ' )(1 ~ ^w*^ 



dz' 



Observe that since ojj is homogeneous of degree 0, for each j 
1, 2, d — 1 and small |z'| we have 



<9z ,c 



< c- 



fc = l,2,.. 



Using this, ((15)) , ((161) , Lemma lLTl and applying Lemma lL9l we obtain 

6 

£ |/ 2 ;>)| < a- 2 J^x), 

mtO,meX d k=l 
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where 



eV 2 <\z'\<C 

£ l/2<| z '|<2 e l/2 

r i i 

A 3 x = — - -dz', 

J \z'\ d \x-z\ a 

£ 1/2<| z '|< 2£ 1/2 
£ l/2<| z /|< 2£ l/2 



and 



f 1 1 
An easy calculation shows that 



d+l 



dz'. 



It 



( £ -3/2 



A k (x)dx < C 



d = 2, 
e _1 |ln£|, d = 3, 
d>4, 



where we used Fubini's theorem to change the volume and surface 
integrations. Using this we obtain 

, £ 1/2 , d = 2, 
£ J |J^(x)|^<C^|ki£|, d = 3, 



d>4. 



Combining together the estimates for L' J and I 2 '^, and using (fl~4|) we 
arrive at 



(17) 



(V /2 , d = 2, 
e|lne|, d = 3, 
e, d > 4. 
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Step 5. U estimates. By virtue of ® we have sup \\u E - WoIIl~(d) < 00 

£>0 

and hence 

(18) ^ \u E (x) - u (x)\dx < Ce, 

D\D f 

which combined with ((17[) gives the claim when p = 1. 

Now for 1 < p < 00 using the boundedness of u £ - u we obtain 

f\u,-u^ X <cj\u l -u < ,m = Q\u,-u Am , 

D D 

hence 

(19) \\u £ - Uolb(D) < Q\u £ - u \\^ p {D y 
Theorem II .21 is proved. □ 



Proof of Theorem 11.41 It is known that 

v £ (x) = u £ (x) - J u £ {y)do{y)- J N(x,y)F £ (y)dy = 

T D 

J N(x,y)g £ (y)do(y). 
r 

By this, the proof of the theorem basically follows from the proof of 
Theorem II .21 by simple modification. Uniform boundedness of v £ (x) 
with respect to x e D and e > now follows from the estimate of 
\N(x, y)\ provided by Lemma fl. 101 Next, instead of ((12)) we consider 
D £ := {x e D : d(x) > e 2 }. It is left to replace P(x,y) by N(x,y), and 
instead of Lemma 11.71 use estimates of Lemma 11.101 in the proof of 
Theorem ll.2l Note that a better convergence rate in comparison with 
the Dirichlet problem is due to lower singularity of the kernel N(x, y) 
than that of P(x, y). □ 



Proof of Theorem 11.61 By the second part of Lemma ll.lOl we see that 
the gradient of the Neumann matrix with respect to x variable, which 
is inside the domain, enjoys almost the same regularity properties as 
that of the Poisson kernel provided by Lemma lL7l But this regularity 
is enough to repeat the steps of the proof of Theorem 11.21 up to Step 
5, and to obtain the same estimates in the region D £ which is away 
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from the boundary. To complete the proof of the theorem for p = 1 
we need to prove the analogue of ((TBI) . Here we use Lemma 11.111 
Keeping the same notations as in the proof of Theorem 11.21 for each 
small t > we have 

£ 

(20) j \V(v £ - v )(x)\dx < C T j J±-<C T j^L = C T e l -\ 

D\D t D\D f 

which proves the case p = 1. Now for 1 < p < oo, take r > p. By the 
Holder's inequality we obtain 

(21) \\V(v e - v Q )\y {D) < \N(v £ - v Q )\\ a L [ (D) \\V(v E - vo)^ 

where l/p = a r + (1 - a r )/r, and a r e [0,1]. From which we conclude 

lr-p 

«r = 7- 

pr-1 

From Lemma fl .111 we have ||V(t7 £ - ^o)!!^^) ^ C T , where C T depends 
on a small parameter in the lemma. The bound for L p -norm now 
follows form the case p = 1 and the fact that lim a r = l/p. □ 

r— >oo 

3. Optimality: proof of Theorem 11.31 

Throughout this section instead of systems we will consider equa- 
tions, so the operator L is considered only in the case N = 1. We 
begin with a simple lemma. 

Lemma 3.1. (Concentration near the boundary) Let u be the solution to 
the Dirichlet problem for the operator L in domain D with boundary data 
g : R rf — > C which is Lipschitz with constant Lip(g). 

Then there exist constants C\, C2 depending on dimension, domain, 
operator, and independent of g, so that for any x £ D, E, e T, and small 
enough 5 > one has 

\u{x)-g^)\<C l bUp{g) + ^\\g\U, 
provided \x - £| < C26. 

Proof. By the Poisson representation we have 

u(x) = J P(x, y)g{y)do{y), xeD. 
r 

Fix£ G Tandx £ D. If < 5/2, and |£-y| > 6 where y £ T, then 
clearly \x - y\ > 5/2. Using this, ((4]), Lemma [TTTI for a = 0, and the fact 



APPLICATIONS OF FOURIER ANALYSIS IN HOMOGENIZATION AND BOUNDARY LAYERUZ 

that the Poisson kernel has integral equal to one over the boundary 
T, we obtain 

(22) \u(x) - g(£)\ = | J \P(x, y)\[g(y) - g(£)]<fo(y)| < 

J \P(x,y)\\g(y)-g(Z)\do(y) + J \P(x,y)\\g(y)-g(Z)\do(y)< 

ly-{|<5 ly-£l>5 

C5Lip(g) + C\\g\\ 00 d(x) J j^j, 

ly-SI>s 

where the constant C is determined by the Poisson kernel. The last 
integral is estimated as in HI and uniformly with respect to x we 
obtain 

diam(D) 



C da(y) r 1 
J \x - y\ d ~ J r d ~ 5' 



y\ a 

ly-4l>6 5/2 

It is left to take x so that |x — £| < C^b, where C2 is a sufficiently small 
constant independent of x G D, £ G T, and g, hence the claim. □ 

The next Lemma is essentially the Weyl's equidistribution theorem, 
in our case concerning equidistribution of the scaled surfaces modulo 
one. 

Definition 3.1. For x, y G R rf we say that they are equal modulo one, and 
write x = y( mod 1) ifx - y G 

If x G R d , by x mod 1 we denote the unique point y in the unit 
torus of R d which is equal to x modulo one. 

Lemma 3.2. (Equidistribution of scaled surfaces) Suppose T is a uniformly 
convex smooth hypersurface in R rf (d > 2). Then for any Riemann integrable 
function g : T d — > IR one has 

(23) J g(x)dx = lim j^jfj J g(fy)My), 

T d r 
where "Kd-i denotes (d - \)-dimensional Hausdorjf measure. 

Proof We first prove the Lemma for smooth functions. Suppose 
g G C 0O (R rf ) and is one periodic. Then 



g(x) = £ c m e lnim - x , x g T d , 
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which converges absolutely. Plugging this expansion into ([23]) we 
see that is is enough to prove that 



a A :=Yj C ™ f e 2mAnty do(y) 



m*0 



converges to 0, as A — > oo. Denote by cr(£) the Fourier transform 
of the surface measure o. The following estimate is known (see ||S]], 
chapter VIII, Theorem 1) 

m\ < cier (rf - i)/2 . 

Using this estimate we obtain 

\Cm\ 



\a A 



< 



^|cJF(A m) |<CA-<«V^ 

m*0 m*0 



The last sum converges, due to smoothness of h, and thus the claim 
follows for smooth functions. 

Now if g is a characteristic function of some ball in the unit torus, 
then it is easy to see that there exist a sequence of smooth functions 
/„ and F n , n = 1, 2, ... so that 

1. f n (x) < g(x) < F n (x), x e R d , 

2. lim f (F n (x) - f n (x))dx = 0, 



n^oa ' 



from which it follows that (|23|) holds true for characteristic functions 
of balls. Clearly it will hold true also for their linear combinations, 
i.e. step-functions. Now observe that when g is Riemann integrable 
function, then the same pointwise bounds from above and below 
hold true by means of step-functions, hence the statement □ 

Applying Lemma [3721 to characteristic functions we obtain the fol- 
lowing result. 

Corollary 3.3. Let T be as above, and A c T d be a ball. Then 

(A) - lim <Kd -^y eT - A y mod 

Now we are ready to complete the proof of Theorem 11.31 
Proof of Theorem 11.31 Without loss of generality we may assume 
that the boundary data g has mean value 0, and hence Uq = 0. 

By the Poisson representation we have 

(24) u £ (x) = J P(x,y)g(y/e)do(y), x e D. 
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If g = 0, then we are done, otherwise denote E := {x G T d : \g(x)\ > 
l/2||glU}. Clearly E is an open set, and by passing to a subset of 
positive measure, we may assume that E is a ball. 

Due to Corollary 13.31 there exists a constant Cq > so that for all 
£' > small enough one has 

^-^K A ^{y e T : \g £ (y)\ > \\\g\U > l/2/i(JE). 
Now fix y G T, so that \g £ (y)\ > 1/211^-Hoo, and apply Lemma [3TT1 with 

l llglU l ejlglU 
8dL/p(g,) 8d Lipig)' 

We obtain 

(25) \u E (x) - g £ (y)\ < HglU/4, if x G D, and |x - y| < Ce, 

where the constant C is independent ofe. Since |g e (y)| > llglU/2 on a 
fixed portion of the boundary for all small enough e > 0, inequality 
(|25|) implies that on a fixed portion of the strip B e := {x G B : d(x) < Ct:} 
one has \u £ (x)\ > ||g|U/4, where x G B E . 

Now for 1 < p < oo taking the 1/ norm of u e only on this strip we 
obtain ||w £ || LP ( B ) > Ce 1/p ||g|| 00 , which proves the theorem. □ 

We remark here that Theorem |1.3l gives sharp bounds on convergence 
rate of the homogenization process in dimensions 4 and higher, and 
nearly sharp in dimension 3. For d = 2, and p = 1 we give an example 
for which the convergence rate is exactly e 1 ^ 2 . 

Example (d=2). Let B be the unit disc of R 2 , and g(xi,x 2 ) = e 2mXl . 
Note that g is one periodic and has mean value in the unit torus. 
Consider the following problem: 

Au e (x) = 0, x G B, 

u e(x) = g(x/e), ^ e dB. 

To estimate u e on B we proceed using the method of stationary 
phase (see e.g. see [S|, chapter VIII). Let P(x, y), where |x| < 1, \y\ = 1 
be the Poisson kernel for the Laplace operator in B. We will consider 
u e (x) only at the points |x| < 1/2 where P(x, y) is a smooth function 
with bounded derivatives. Observe that the only critical points of the 
phase function g are north and south poles of the disc, i.e. n+ := (0, 1) 
and n_ := (0, -1). It is also clear that these are non-degenerate critical 
points. Hence we can invoke the principle of stationary phase (see 
ISl , chapter VIII, Prop. 6) and obtain 

u £ (x) = Ce y2 [P(x,n + )e^ +P(x,n_)e"^] + 0(e 3/2 ), 
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where 0(£ 3 ^ 2 ) is uniform with respect to \x\ < 1 /2. Now to see that 
the two terms in the parentheses do not cancel, it is enough to restrict 
x to {x = (xi,x 2 ) € B : \x\ < 1/2, 1/4 < x 2 < 1/2}. Considering 
u e on this subset we see that ||w £ || L i (B) > Ce 1/2 , which proves that the 
convergence rate provided by Theorem ll.3l in the case p = 1 and d = 2 
is optimal. 
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